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Abstract 



The gravitational collapse of a thick cylindrical shell of dust matter is investi- 
gated. It is found that a spacetime singularity forms on the symmetry axis and that 
it is necessarily naked, i.e., observable in principle. We propose a physically reason- 
able boundary condition at this naked singularity to construct the solution including 
its causal future. This boundary condition enables us to construct the unique con- 
tinuation of spacetime beyond the naked singularity and ensures that the dust shell 
passes through the naked singularity. When the cylindrical shell leaves its symmetry 
axis away, the naked singularity disappears, and regularity is recovered. We construct 
numerical solutions with this feature. This result implies that the gravity produced 
by a thick cylindrical shell of dust is too weak to bind the shell even if it engenders 
the formation of a curvature singularity which is so strong as to satisfy the limiting 
focusing condition. For this reason, this naked singularity is very weak in the extended 
spacetime; the metric tensor is C 1 " even at the naked singularity, and the extended 
spacetime is complete for almost all geodesies. This feature is also seen for singular 
hypersurfaces. Such an extended spacetime can be regarded as phenomenological in 
the sense that it is valid if the relevant microphysics length scale is sufficiently small 
compared to the scale of interest. 
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§1. Introduction 



General relativity predicts the gravitational collapse of very massive objects and the for- 
mation of spacetime singularities which will be accompanied by the divergence of physical 
quantities, such as the energy density, stress and spacetime curvature. 1 - ) An important is- 
sue related to the attribute of spacetime singularities is known as the cosmic censorship 
hypothesis. 2 ^ Roughly speaking, this hypothesis asserts the impossibility of the formation 
of an observable spacetime singularity, called a naked singularity in our universe. How- 
ever, preceding studies have revealed that there are many possible counterexamples to this 
hypothesis. 3 ^ 13 ) 

Although the appearance of spacetime singularities is an important consequence of gen- 
eral relativity, spacetime singularities themselves cannot be described in a physically reliable 
manner within the framework of general relativity, because general relativity becomes in- 
valid beyond the quantum gravity scale. Even in the neighborhood of spacetime singularities, 
quantum effects in gravity should be important and in this sense, the neighborhood of sin- 
gularities is also regarded as beyond the applicability of general relativity. In light of this 
observation, two of the present authors introduced the concept of a spacetime border. 15 ^ 
Roughly speaking, a spacetime border is the spacetime region in which general relativity is 
not applicable due to the high curvature of spacetime within it. The implication of a possible 
counterexample to the cosmic censorship hypothesis would be that we can directly observe 
physical processes at the spacetime borders in our universe, other than the Big Bang. 

Some of the known analytic solutions of the Einstein equations for naked singular space- 
times have been extensively investigated, and the structures and strengths of the correspond- 
ing naked singularities have been revealed. 3 ) -9 )' 17 ^ 24 )' 34 ) To determine whether a naked sin- 
gularity is null or timelike and how a naked singularity is observed (such as in gravitational 
redshift), we have to take the limit to a naked singularity from its causal future. This is an 
inevitable dilemma. We want to predict the observational signature of naked singularities, 
while naked singularities are what prevent us from constructing the unique time development 
of the spacetime beyond them. In order to construct a solution of the Einstein equations 
for the causal future of a naked singularity, the boundary condition at the naked singularity 
must be specified, because spacetime singularities are boundaries of the spacetime. 

It might be thought of that a unique boundary condition at naked singularities cannot be 
selected from a large, maybe even infinite number of possible conditions, because we do not 
understand the laws of physics which describe physical processes in the spacetime border 
covering a naked singularity. However, we should note that for all of the known naked 
solutions, some boundary conditions at naked singularities have been implicitly imposed 
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without understanding the physical laws applicable to spacetime borders. One method for 
obtaining solutions for the causal future of naked singularities formed from regular initial 
data is to assume the extension across the Cauchy horizon associated with them to be 
as smooth as possible. A typical example is the subcritical Reissner-Nordstrom solution. 
Many examples of naked singular spacetimes are in this category, in which analyticity on 
the Cauchy horizon is required. In other words, many known naked singular solutions are 
obtained by imposing boundary conditions which make possible smoothness on its Cauchy 
horizon, although nobody yet knows whether this is valid from the point of view of a quantum 
theory of gravity. It is important to elucidate the principle regarding the determination of 
boundary conditions at naked singularities in known analytic solutions, but this is beyond 
the scope of this paper. 

An exceptional example of naked singularities is the class of singular hypersurfaces on 
which some scalar polynomials of the curvature tensor diverge. It is believed that they can 
be described in a reliable manner with the junction condition. 14 ) In fact, it seems that many 
people use a singular hypersurface without realizing that this introduces a curvature singu- 
larity. It is not difficult to understand why this is possible. If the hypersurface is actually 
infinitesimally thin, of course the energy density diverges, and therefore it corresponds to 
a curvature singularity. However, to derive the junction condition, we first assume a finite 
thickness e, integrate the Einstein equation along the normal direction, and then take the 
limit e — > 0, keeping the surface density and pressure constant. After taking this limit, 
the metric tensor can still be locally C l ~, in which case geodesies are uniquely extended 
across the singular hypersurface. 1 -* A singular hypersurface is regarding as describing the 
gravitational effects of a hypersurface of very small but finite thickness. Hence, a singular 
hypersurface can be regarded as providing a phenomenological description of a physical hy- 
persurface of small but finite thickness. This is also the reason why the singular hypersurface 
formalism should be capable of treating gravitational effects of a hypersuface whose thick- 
ness is as small as the quantum gravity scale, such as a D-brane. 25 ) To proceed with this 
approach, if we can find a C 1 " extension of spacetime beyond a naked singularity, it could be 
regarded as a phenomenological description of something whose length scale is much smaller 
than the scale of interest, because of the unique extendibility of geodesies across singularities. 

In this paper, we investigate the gravitational collapse of a cylindrically symmetric thick 
dust shell surrounding a vacuum interior, in other words, a hollow cylinder composed of 
dust matter. Such a configuration of matter causes the formation of a string-like naked 
singularity of codimension two along the symmetry axis. 26 )' 27 ) We propose a physically 
reasonable boundary condition at this naked singularity. It is difficult to analytically obtain 
solutions for this system. Therefore we numerically solve the basic equations. The obtained 
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solution is a C extension of spacetime, and geodesies uniquely extend across the naked 
singularity, except for that which is on the string-like naked singularity. This system is not 
asymptotically flat, and thus it alone does not constitute a counterexample to weak cosmic 
censorship. However, this model might describe a portion of the spindle collapse studied by 
Shapiro and Teukolsky, which is an interesting possibility as a counterexample to the cosmic 
censorship hypothesis. 13 - 1 

This paper is organized as follows. In §2, we give the basic equations for the gravitational 
collapse of dust matter in a spacetime with whole cylinder symmetry. In §3, we propose 
physically reasonable boundary conditions on metric and matter variables on the symmetry 
axis where a naked singularity is caused to form by the cylindrical gravitational collapse of 
a thick dust shell. Then, we discuss the physical meaning of such boundary conditions. In 
§4, we present a numerical solution that describes the gravitational collapse of a thick dust 
shell. In §5, we show that almost all of the causal geodesies are complete in the spacetime 
considered in this paper. In §6, the strength of the curvature at the naked singularity is 
investigated. Finally, §7 is devoted to summary and discussion. In Appendices A and B, we 
prove two facts which are very important in in the investigation carried out in this paper. In 
Appendix C, the Christoffel symbols of the spacetime with whole cylinder symmetry in the 
Cartesian coordinate system are given, and in Appendix D, the components of the Riemann 
tensor in the cylindrical coordinate system are given. 

In this paper, we adopt units in which c = 1. Greek indices denote components with 
respect to the coordinate basis. 

§2. Basic equations for a cylindrical dust system 

We consider the spacetime with whole cylinder symmetry 28 )' 29 ) defined by the metric 

ds 2 = e 2 ^ (-dt 2 + dr 2 ) + r 2 f3 2 e" 2 ^d V 2 + e 2 ^dz 2 , (2-1) 

where < r < +oo, < ip < 2-n and — oo < z < +oo constitute the cylindrical coordinate 
system, and 7, ip and (5 are functions of the time coordinate t and the radial coordinate r. 
The coordinate basis vectors d/dep and d/dz are Killing vectors. The coordinates t, r and z 
are all normalized so as to be dimensionless. 

We consider a two-component model of dust matter and assume that one component 
is collapsing while the other component is expanding and that these components do not 
interact with each other. The stress-energy tensors of the collapsing dust, T? v ', and the 
expanding dust, T£ u , are, respectively, written in the forms 

Tr = pu»u u , T» v = fiv^v u , (2-2) 
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where p and p are the rest mass densities, whereas and are the 4-velocities of the 
collapsing and expanding components of the dust matter, respectively. We assume that p 
and p are non-negative so that all of the energy conditions are satisfied. By virtue of the 
assumed symmetry of the spacetime, the 4-velocities and v M can be written in the forms 



;i,tt,o,o) 



[l,v,0,0) 



(2-3) 



where u (< 0) and v (> 0) correspond to the 3-velocities of the collapsing and expanding 
components of the dust matter, respectively. 

In place of the rest mass densities p and p, we introduce the conserved rest mass densities 
D and E defined by 

r/3e^p 



D :-. 
E :-- 



-gpu 



-gpv 



u 



vT 

r(3e^p 

VT^T 2 "' 



(2-4) 
(2-5) 



where g is the determinant of the metric tensor. Then the equations of motion for each 
component of the dust matter, V U T£ U = and V v T? l/ = 0, become 



dt 



D + (uD)' = 0, 
E + (vE)' = 0, 

(1 - u 2 ) u (Jj - jj + ip' - j' 
(1 - v 2 ) \v (ij) - 7 ] + v' - i 



du 

— = u + uu 
dt 

dv 



V + vv 



(2-6) 
(2-7) 

(2-8) 
(2-9) 



where the dot represents the time derivative, while the dash represents the radial derivative. 
It should be noted that the first and second equations represent conservation of the rest 
mass, and for this reason we refer to D and E as the conserved rest mass densities, whereas 
Eqs. (12-81) and (12-91) are the geodesic equations for the constituent particles of the dust 
matter. 

The Einstein equations take the forms 



1 >=\P + 2rf3f3> + r 2 ((3' 2 -f3 2 ) 



r(3 {(3 + r(3') ( ^ + tp' 2 ) - 2r 2 (3(5^' 



+ 2pp' + r ( 2p' 2 + [3(3" -0 2 )+r 2 [ (3' (3" - (3(3' 



+ 



VT 



(3 + r(3' + r(3u) + 



8nGe^E 



VT 



(3 + r(3' + r(3v 



(2-10) 
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2 + 2r(3(3' + r 2 ((3' 2 - (3 2 ) 



r 2 (3(3 U 2 + ij' 2 ) - 2r(3 ((3 + r(3') ^ 
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-PP + r (pp - pp') + r 2 (j3P" - P'P') 



P-P"- -P' = —e^ ( DVT^ + EVl^A , 
r r \ J 

7-7 = -0 - ip , 




(2-11) 



(2-12) 



(2-13) 




x ( DVl-u 2 + EVl-v 2 



(2-14) 



Equations (12-101) and (12-111) are the constraint equations, and Eqs. (12- 121) — (12- 141) are the 
evolution equations for the metric variables P, 7 and ijj. 



3.1. Spacetime regularity 

First, we stress that it is possible to construct C 2 ~ solutions for the metric variables, P, 7 
and ifj, and everywhere finite solutions for the matter variables, D, E, u and v, as functions 
of t and r, even if the spacetime singularity is caused to form at r = by the gravitational 
collapse of the cylindrical dust matter. This may seem strange. However, it should be noted 
that it is not a sufficient condition for the realization of regular spacetimes that the metric 
variables be C 2 ~ functions of t and r. As shown below, in order for the spacetime to be 
regular, the first-order radial derivative of the metric variables P', 7' and if)' should vanish 
at r = 0, whereas the matter variables D, E, u and v should vanish at r = 0. In order 
to understand the meaning of these conditions, we introduce a Cartesian coordinate system 
defined by 



with the remaining coordinates, t and z, unchanged. In the case of the regular spacetime, 
the coordinate singularity at r = in the cylindrical coordinate system does not exist in this 
Cartesian coordinate system. 

We have the following components of the metric tensor in this new coordinate basis: 



3. Boundary condition 



x = r cos ip, 



y = r sin <p, 



(3-1) 




(3-2) 



(3-3) 



(3-4) 



7 



9yy 

9zz 



(3 2 e~^ 

2ip 



P 2 



1 sin z tp + 1 



(3-5) 
(3-6) 



with all other components vanishing. It can be easily seen from the above equations that 
(3 2 should be equal to e 27 at r = in order for the components of the metric tensor to be 
single-valued at r = 0. (Note that any value can be assigned to (p at r = 0.) The first and 
second order derivatives of the components of the metric tensor in the Cartesian coordinate 
system with respect to x and y are written in the form 



dx 
dg^u 
dy 

d 2 g„ v 
dx 2 



g 



, 1 9g^ ■ 

cos ip H sm (p, 

r dip 

, . 1 dg^ 
V sinp--— cosp, 

l -gl v {\ + cos 2<p) + ± 



(3-7) 
(3-8) 



gLi 1 - cos2^) - 



9g'^ . ' 

sm lu) 
dip 



1 

+ -T 



dg 



d 2 g» v 
dxdy 



d 2 g^u 
dy 2 



1 f . . _ dg' 



(3-9) 



—a" sin2cp — 
2 ^ Y 2r 



gL sin 2 v - -^i 1 + cos 2( p) 



dip 



cos 2</? + - v, sm 2^ 



dip 



2 <V 



(3-10) 



cos2y?) H 

2r 



5^(1 



cos2y?) + 



sm 2</? 

ay? 



//// 



sin 2^9 



-^(l + cos2^) 



(3-11] 



We can see from the above equations and Eqs. fl3-2l) - fl3-6p that in order for the components 
of the metric tensor g^ in the Cartesian coordinate system to be C 2 ~ functions of x and y, 
the following regularity conditions must be satisfied: 

P 2 7 _ ^2 



lim — 

lim — 

lim 0" 



[finite] , 



[finite] , 
[finite] , 



lim — 

r— >0 r 

lim 7" 



[finite] , 
[finite] , 



lim — 

r— >0 7* 

lim ■0" 



[finite] , 
[finite] . 



(3-12) 

(3-13) 
(3-14) 



Note that the above regularity conditions do not guarantee the continuity of the second 
order derivatives of the metric components at the symmetry axis, r = 0. Although this does 
not necessarily mean that the metric tensor is not C 2 , it might imply that the Cartesian 
coordinate system (13- ip itself is C 2 ~ . However, this Cartesian coordinate system is sufficient 
for our present purpose. 
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We can see from Eqs. (!2-4p and ( 12- 5ft that .D and E vanish on the symmetry axis, r = 0, if 
the rest mass densities p and \i and the metric variables are everywhere finite, and, further, 
if the absolute values of u and v are smaller than unity. The components of the 4-velocities 

and v M in the Cartesian coordinate system are given by 

= — = (l,u cosy, u sin <p, 0) , (3-15) 



f M = — ^= (1, f cosy, f sin y?, 0) . (3-16) 



From the above equations, it is seen that in order for the components of w M and v M in the 
Cartesian coordinate system to be everywhere single- valued, u and v must vanish at r = 0. 
Thus, the regularity condition on the matter variables are 

limD = limi? = and limw = limv = 0. (3-17) 

r^0 r->0 r^0 

3.2. Functional regularity and boundary conditions on metric components 

If a spacetime singularity is formed by the collapse of the cylindrical dust matter, the 
regularity conditions on the metric variables, (I3-I2p and (13-131) . and on the matter variables, 
(13 -171) . will not be satisfied. From Eq. (12-121) . we have 

/3' = -47rGe 7 ^ [d^I^ 2 + EVT^ 2 ) + r [f3 - /3") . (3-18) 

Because we construct C 2 ~ solutions for the metric variables and everywhere finite solutions 
for the matter variables, the above equation gives the following Neumann boundary condition 
on (3 at r = 0: 

fj'\ r=0 = -AnGe 1 ^ (dVY^? + EVT^^j ^ . (3-19) 

With the same procedure, we obtain the following Neumann boundary conditions on ip and 
7 at r = from Eqs. fl2~T0|) and fl2~T4j) : 

SnGe^ / Du 2 Ev 2 \ , . 

= — iTO + 7T=^J „„ ' (3 ' 20) 

^'l^^-^e^ (dVT^ + eVT^ 2 ) . (3-21) 

fj V / r=0 

Here we have used Eq. (13- 19j) to derive Eq. (!3-20p . 

Because it is difficult to analytically construct solutions for cylindrically symmetric space- 
times with dust matter, numerical simulations are necessary to study their detailed behavior. 
The boundary conditions (13- 19j) ( |3-21j) are sufficient to construct numerical solutions for the 
evolution equations (12- 12j) (12- 14j) without any ambiguities, if the matter variables D, E, u 
and v are everywhere finite. 
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There will be three kinds of spacetime singularities formed through the collapse of a 
cylindrical dust matter. The first kind is the so-called shell crossing singularity, whose 
appearance is caused by a caustic of dust matter at some position not on the symmetry axis 
and is necessarily accompanied by the divergence of one or both of the conserved densities 
D and E. However, it is known that shell crossing singularities are so weak that they are 
not physically so significant. For this reason, we do not consider shell crossing singularities 
here. The second kind of spacetime singularity forms on the symmetry axis, r = 0, without 
a caustic of the dust matter. In this case, the matter variables D, E, u and v are everywhere 
finite, even if a spacetime singularity does exist. This is the case of interest in this paper. 
Since this kind of spacetime singularity is formed by the dust matter focused on r = 0, it 
is often called a shell focusing singularity. The third kind of spacetime singularity is also a 
shell focusing singularity, but with a "caustic" of dust matter at r = 0. In this case, D or 
E diverges at r = 0, and thus it is impossible for the metric variables to remain C 2 ~ and 
the matter variables to remain everywhere finite, even if we impose the boundary conditions 
(I3-19I) -( 13-21I) . In this paper, we do not consider this kind of spacetime singularity in detail, 
but we briefly discuss it in §7. If (3 vanishes at some point, it would seem that another 
kind of spacetime singularity is formed there, as several terms in the Einstein equations 
(12-10p - (l2-14p are proportional to (3^ 1 . However, as shown in Appendix A, this possibility is 
excluded by the boundary conditions on the matter variables. 

Hereafter, we focus on the second kind of spacetime singularity. Now we can see how a 
spacetime singularity forms on the symmetry axis r = through the gravitational collapse of 
the dust matter while the metric variables (3, 7 and ip remain C 2 ~ and the matter variables 
D, E, u and v remain everywhere finite. Since we are interested in gravitational collapse, 
we consider the initial data with the only collapsing dust (D,u). To avoid the appearance 
of a caustic before the formation of the second kind of spacetime singularity, we assume 
that the 3-velocity field u is initially a monotonically increasing function of r. (Here note 
that we have u < 0, by definition.) If p initially does not vanish on the symmetry axis 
r = with such an initial 3-velocity field, the regularity condition on the 3-velocity field 
u in (13-171) cannot be satisfied. Thus, we have to assume that the initial configuration of 
dust matter is a cylindrical thick shell, i.e., a hollow cylinder. (The case of a solid cylinder 
is discussed in §7.) When a collapsing thick dust shell reaches the symmetry axis, r = 0, 
D becomes non- vanishing at r = 0. Then the rest mass density p of the collapsing dust 
becomes infinite, because p is proportional to D/r, as seen from Eq. f l2-4j) . Because the Ricci 
scalar is proportional to p + p as given by the Einstein equations, the divergence of p implies 
the formation of a scalar polynomial (s.p.) curvature singularity. ^ It is also easily seen from 
Eqs. (13- 19j) (I3-21|) that if D does not vanish at r = 0, then /?', 7' and ip' do not vanish at 
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r = 0, or equivalently, the regularity of spacetime is broken. 

3.3. Boundary conditions on dust matter 

It should be noted that we still have freedom in the choice of the boundary conditions 
on the matter variables D, E, u and v imposed at r = 0. Here we impose the following 
boundary conditions on the matter variables at the symmetry axis, r = 0: 

E(t, 0) = D(t, 0) and v{t, 0) = -u{t, 0). (3-22) 

The above boundary conditions physically mean that once the collapsing dust reaches the 
symmetry axis r = 0, the same amount of the expanding dust is shot from r = with the 
same speed as the collapsing dust matter. 

In order to impose the above boundary conditions in numerical simulations, it is conve- 
nient to extend the domain < r < oo to — oo < r < oo; we call the original domain the 
physical domain and the additional domain, — oo < r < 0, the fictitious domain. We de- 
fine the metric variables in the fictitious domain such that these variables possess reflection 
symmetry with respect to r = 0, i.e., 

7 (t, - r ) = 7(t, r) and ip(t, -r) = if>(t, r). (3-23) 

Then we solve Eqs. (12-61) and (I2-8|) for the variables of the collapsing dust (D,u) in the 
fictitious domain, r < 0, as well as in the physical domain, r > 0. For the variables of the 
expanding dust (E,v), we impose the conditions 

E(t, r) = D(t, -r) and v(t, r) = -u(t, -r). (3-24) 

The variables E and v determined by the above condition automatically satisfy Eqs. (I2-7|) 
and (J2U). 

Here it should be noted that in the extended domain, — oo < r < oo, the metric variables 
are no longer C 2 ~ , even though they are regular with respect to t and r in the physical 
domain, < r < oo. If a spacetime singularity forms, 7 ' and ip' do not vanish at the 
symmetry axis, r = 0. Then the reflection symmetry condition (13-231) leads discontinuities 
of 7' and ip' at r = in the extended domain. This low differentiability in the extended 
domain causes a decrease in the accuracy of numerical simulations, and thus we need a 
careful treatment to maintain good numerical accuracy. 

3.4. Physical consequences of the boundary conditions 

The above conditions on the matter variables lead the following three physically favorable 
consequences. 
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(i) C l metric 

It should be noted that the boundary condition (I3-24|) guarantees that the metric tensor 
is defined even at the spacetime singularity at r = 0. From Eq. (12- lip , we find the relation 

£ f£) = -arfW^V f _£^ + _£U atr = (3-25) 
dt\i3j \P J \^T=tf s/Y^fl) ( ' 

is satisfied. Before the singularity formation, e 27 / (3 2 should be unity. (Hereafter, without loss 
of generality, we assume e 1 / (3 = 1 at r = before the singularity formation.) Otherwise, the 
symmetry axis, r = 0, would be conically singular even before the thick dust shell collapses 
to the symmetry axis, r = 0. Thus the condition (13-241) and the above equation always 
guarantee the relation 

f_ = 1 at r = 0. (3-26) 

As previously shown, this implies that the metric tensor is finite and single-valued, even 
at the spacetime singularity, or, in other words, it is at least C°. Further, it is not such a 
difficult task to show that the components of the metric tensor in the Cartesian coordinate 
system are C 1_ , i.e., locally Lipschitz functions. 

(ii) Completeness of radial geodesies 

As shown in §5, the radial geodesies that intersect the spacetime singularity at r = 
from the domain r > can be extended across the spacetime singularity. This implies that 
the trajectories of dust particles can be extended across the spacetime singularity that they 
cause to form. 

(iii) Conservation of the total rest mass in the physical domain 

From Eq. ()2-6p . we find that the rest mass of collapsing dust per unit Killing length a 
along ^-coordinate in the extended domain, — oo < r < oo, is constant, where a is given by 

/oo 
D{t,r)dr. (3-27) 
-oo 

By the condition (13-241) . this can be rewritten in the form 

/■oo 

a= [D(t,r) + E(t,r)]dr. (3-28) 
Jo 

The above equation implies the conservation of the total rest mass in the physical domain, 
r > 0. 

As demonstrated by our numerical results presented in the next section, if we set the 
initial conditions as D > and u < with E — = v in the physical domain, all of the 
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Fictitious Domain 




Physical Domain 



dust matter 



t matter 



► r 



Initial hypersurface 



Fig. 1. A schematic diagram of the gravitational collapse of a thick cylindrical shell of dust. The 
collapsing dust forms a naked s.p. singularity, and then dust particles pass through the naked 
singularity created by the dust itself. The regularity at the symmetry axis, r = 0, is recovered 
after the dust leaves the naked singularity. 

dust matter will enter into the fictitious domain, r < 0, from the physical domain, r > 
(see Fig. 1). By the condition f)3-24p . the collapsing dust (D,u) in the fictitious domain is 
equivalent to the expanding dust (E,v) in the physical domain. Further, it should be noted 
that the completeness of radial geodesies implies that the dust matter can be regarded as 
a cold gas composed of collisionless particles; the collapsing dust (D, u) passes through the 
spacetime singularity that it creates and then leaves the symmetry axis r = as expanding 
dust (E,v). In this sense, it is natural to call the condition fj3-24j) the passing-through 
condition. As mentioned above, when the dust matter passes through the symmetry axis, 
the rest mass densities of dust, p and p and, thus, the Ricci scalar diverge there. 



§4. Numerical simulation 



4.1. 



Initial conditions 



We choose the initial conserved rest mass densities and 3-velocity fields to be 



D 



15a 



[r - 1(1 - w)} 2 [r - 1(1 + w)f , 



(4-1) 
(4-2) 



32t™ 5 / 5 



E = 0, 
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« = - A /l-exp(~), (4-3) 
v = 0, (4-4) 

for /(l — w) < r < Z(l + w) and vanishing elsewhere, where a is defined by Eq. (13-281) . I 
and w (< 1) are positive parameters controling the location and thickness of the dust shell, 
respectively, and v is a parameter controling the gradient of the 3-velocity field, u. 

We determine the initial data of the metric variables /3, 7 and ip and their time derivatives 
in the following manner. First, we set (5 = 1 and (3 — 0. Then the constraint equations 
(12401) and (EED) become 

1 =r%l) + — , (4-5) 



v 7 ! 



^-^r^- < 4 ' 6 > 

We have to integrate Eq. (14-51) numerically in order to obtain 7, while Eq. (14-61) directly 
gives its time derivative. We are interested in an initial situation that is as close to the static 
configuration as possible with a non-vanishing initial collapsing velocity. Therefore we set 
■0 = 0. In order to determine the initial data of ip, we use Eq. (12-141) with ip = 0, i.e., 

0" = -- (V + ATrGe^DVT^^j . (4-7) 

We numerically integrated Eqs. (14-51) and (14- 7p simultaneously outward from r = by 
imposing the boundary conditions 7| r =o = and V|r=o = = i/j'\ r =o, since we assume 
that there is no spacetime singularity in the initial data. 

The vacuum region in the initial data obtained with the procedure described above agrees 
with the Levi-Civita solution, which is the unique solution for a vacuum static spacetime 
with whole cylinder symmetry, 

P — 1, 7 = /t 2 lnr + A and ip — — ftlnr, (4-8) 

where k and A are parameters that characterize this solution. Integrating Eq. (I4-7|) . we find 
that k vanishes for r < 1(1 — w), whereas for r > 1(1 + w), we have 

k = 47tG / dre^'DVl - u 2 . (4-9) 

Because D is non-negative, k is positive in the domain r > 1(1 +w). 
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Fig. 2. Several snapshots of (3 in the physical domain, r > 0. The horizontal axis represents the 
radial coordinate, r. 
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4.2. Evolution 

In order to solve the evolution equations for the metric variables (I2-I3j) — f l2-14j) numerically, 
we replaced the spatial derivatives by finite differences and adopted the MacCormack scheme 
for the time integration. In order to solve the evolution equations for the dust matter, 
we adopted the method proposed by Shapiro and Teukolsky, 13 ) with which we follow the 
motion of constituent mass shells of the dust matter moving along timelike geodesies. For 
the integration of the geodesic equations, we also adopted the MacCormack scheme. Then 
we constructed the conserved rest mass densities D and E and the velocity fields u and v 












: 1 \J 

t=0 




t=4 




/ s 




(=12 




(=16 



5 10 15 20 



r 

g. 3. The same as Fig. 2, but for 7. 
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Fig. 4. The same as Fig. 2, but for 7. 
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Fig. 5. Several snapshots of the conserved rest mass density D in the extended domain, — 00 < 
r < 00. The horizontal axis represents the radial coordinate, r. 

from the positions and velocities of these mass shells. We assumed that all of the mass shells 
have the same rest mass per unit Killing length along z-coordinate. A linear combination of 
the constraint equations ( I2T0I) and ( l2Tlj) was used to check the accuracy. 

We display examples of numerical solutions for the metric variables /3 — 1, 7 and ip in 
the physical domain and the matter variables D and u in the extended domain in Figs. 2-6, 
respectively, where we have chosen the parameter values I — 1, w — 0.5, a = 10~ 5 and 
v = 10~ 2 . The numerically covered domain is —20 < r < 20. The numbers of spatial grid 
points and mass shells are 5 x 10 3 and 250, respectively. In this numerical simulation, the 
spatially averaged relative error estimated from the constraint equations was always less 
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Fig. 6. The same as Fig. 5, but for the velocity field u. 



than 10~ 8 . Another criterion for the numerical accuracy is the difference of /3e _7 | r= o from 
unity. This difference was also always of order 10~ 8 in this simulation. Because \/3 — 1| and 
7 are at most of order 10~ 5 , the order of the relative error is understood to be 10 -3 . 

The boundary condition at the outer numerical boundary, r = 20, is that none of the 
variables change there. This simplistic boundary condition is valid because we are interested 
in the evolution of the central region, and we stop the numerical simulations before the 
influence of the non-trivial initial configurations of the metric and matter variables reach the 
numerical boundary, r = 20. 

It is seen from Figs. 2-4 that the metric variables f3, 7 and if) are continuous and smooth 
in the numerically covered region. Further, it is seen from Figs. 5 and 6 that the conserved 
rest mass density D is everywhere finite and that the absolute values of the velocity field 
u is always smaller than unity. However, as mentioned above, when the dust shell reaches 
the symmetry axis, r = 0, the rest mass density p diverges there, and thus a s.p. curvature 
singularity 1 -* forms on the symmetry axis, r = 0. The dust matter reaches the symmetry 
axis at t ~ 3.0 and forms a spacetime singularity. Then, all of the dust matter leaves the 
spacetime singularity, and regularity at the symmetry axis is recovered at t ~ 12.3 (see 
Fig. 5). The propagations of f3 and 7 represent the gravitational radiation generated by the 
gravitational collapse of the dust matter. 

Because even after the appearance of the spacetime singularity, (3, 7 and if> are everywhere 
finite and continuous, there is a causal future of this spacetime singularity, and there exists 
a Cauchy horizon associated with it. We can thus conclude that a naked singularity has 
indeed formed in the spacetime constructed by this numerical simulation. 
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}5. Completeness of causal geodesies 



In the spacetime of interest, since the metric tensor is C 1 " even in the neighborhood of 
the spacetime singularity, the parallel relation between any pair of vectors can be defined 
everywhere. A pair of spacelike vectors or timelike vectors p^ 1 and g M are parallel if and only 
if the equality 

|gg£ = ±1 (5-1) 

is satisfied, where we have +1 for the pair of the spacelike vectors and —1 for the timelike 
vectors. In the case of a pair of null vectors and Z M , they are parallel if and only if 
g^ish^l" — and the product of their time components are positive. Further, in the spacetime 
of interest, the Christoffel symbols in the Cartesian coordinate system defined by Eq. (13- ip 
have limits with directional dependence for r — > 0, but these limits are finite (see Appendix 
C). By virtue of these two properties of the spacetime singularity at r = 0, almost all causal 
geodesies are complete, even if they intersect this spacetime singularity. We show how to 
construct complete solutions for the geodesic equations below. 

In the coordinate system (12- 1[) . the geodesic equations take the forms 



d 

dr 



,2( 7 -V0 



|] = (7-^)e 



,2( 7 -V0 



fdtV fdr\ 



dz 



(5-2) 
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d 
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2 fdr\ 2 ' 
dr J \dr J 



dtV 



r 



\dr J 



dr J 

dz\ 
Tfr) 



0. 



dr 



(5-3) 
(5-4) 
(5-5) 



where r is the affine parameter. The third and fourth equations can be easily integrated 
once, and we have 



dr 



e - " — = P, 
dr 



(5-6) 



where L and P are constants of integration. There are three categories of causal geodesies, 
which behave in very different manners near the spacetime singularity. The first category 
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consists of the causal geodesies with r = = dr/dr. We call these geodesies central geodesies. 
The second category consists of the causal geodesies with vanishing L but non-vanishing 
dr/dr. We call causal geodesies in this second category radial geodesies. The third category 
consists of geodesies which have non-vanishing L. We call these non-radial geodesies. As 
shown in Appendix B, non-radial geodesies cannot reach the symmetry axis, r = 0. For this 
reason, hereafter we focus on the central and radial geodesies. 

5.1. Radial geodesies 

First, we consider radial geodesies. Hereafter we assume that all causal geodesies are 
future directed, i.e., their time components are positive. An extension of a radial geodesic 
beyond the spacetime singularity at r = is also a radial geodesic. Thus we consider 
two radial geodesies, denoted by fe_ and k + . We assume that fc_ intersects the spacetime 
singularity at t = to and that k + leaves the spacetime singularity at t = to Thus k + U k_ 
consists of a continuous curve. In order for k + to be a unique extension of k_ beyond the 
spacetime singularity, the tangent vectors of these two geodesies should be identical at the 
spacetime singularity, since the geodesic tangent vector is parallelly transported along the 
geodesic, by definition. If this condition uniquely determines k + for any we can conclude 
that all of the radial geodesies are complete in this spacetime. 

In general, a radial geodesic is specified by = [constant]. We denote the angle coordi- 
nates of k± by tp±. The components of their tangent vectors, k±, in the cylindrical coordinate 
system are written in the form 

fc£ = (yJ-eW-r) X + V£ + e-*iPL V ± , 0, e^P ± ) , (5-7) 

V / tripz 

where V± are functions of the affine parameter r, P± is constant, and \ is —1 if k± are 
timelike geodesies, while x vanishes if k± are null. Here we adjust the affine parameter r 
so that the geodesic fc_ intersects the spacetime singularity at r = 0, while k + leaves the 
spacetime singularity at r = 0, or, in other words, a negative r is assigned to fc_ and a 
positive r is assigned to k + . By assumption, we have 

lim V_ < and lim V + > 0. (5-8) 

In order to see the continuity of the tangent vectors at the spacetime singularity, the Carte- 
sian coordinate system defined by Eq. (13-11) is useful, because the Cartesian coordinate 
system covers the symmetry axis, r = 0, whereas the cylindrical coordinate system does not, 
and the components of vectors in this coordinate system cannot be specified at r = 0. The 
components of the tangent vectors of k± in the Cartesian coordinate system are therefore 
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given by 



k^ = (yJ-e 2 W-^x + VI + e- 2 ^Pl V±coscp±, 14sin^±, e~ 2lp P ± j 

\ / txyz 



(5-9) 



Since the metric variables (5, 7 and ip are C 2 ~ functions of t and r, the quantities V± have 
finite limits as r — * 0. The components of the tangent vectors k± approach identical values 
in the limit that the spacetime singularity is approached, r — > 0, if and only if the following 
conditions are satisfied of these: 



P + = P. 



lim V+ 

T->0 



lim V_ 



and ip + = f - + 7T. 



(5-10) 



The first of these equations is obtained from the continuity of the z-component, whereas 
the second one is obtained from the continuity of the t-component and Eq. (15-8p . The final 
one is derived from the continuities of the x- and y-components. Here it should be noted 
that because the metric tensor is defined even at the spacetime singularity at r = 0, the 
fact that the limiting values of the components k± obtained as r — > are identical implies 
that the tangent vectors themselves coincide at the spacetime singularity. Since the first 
order derivatives of 7 and ip are C 1_ functions of t and r in the domain < r < 00, the 
standard existence theorems for ordinary differential equations can be applied to the geodesic 
equations for radial geodesies k + , 



dV+ 

dr 
dt 

dr 
dr 



3 2W-7) 



(2^ - i)e~^Pl + (V> - 7)X 



+ 2^-7), 



+ • 



(5-11) 
(5-12) 
(5-13) 



and Eq. (!5-6[) with L = and P = P + , in order to show that the geodesic k + is uniquely 
determined by P + , V + and y? + at (t, r) = (t , 0). Therefore we conclude that k + U fc_ is the 
unique extension of the geodesic fc_ beyond the spacetime singularity. 

5.2. Central geodesies 

In order to elucidate the extendibility of the central geodesies, first we consider a causal 
curve along r = 0. By definition, the tangent vector to this curve is given by 



Thus we have 



dx 11 dx v 
tJ,v ~d\~d\ 



dx^ 
~d\ 



-, 0, 0, - 
dX' dX 



dt\ 
dXJ 



ip'e 



-27+41/) 



dz V 



COS(y9, 



(5-14) 



(5-15) 
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dx^ dx 
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dX 
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dz 
dX 



) 



2 



^ dr dX 



(Y - 



sin v 5 , 



(5-16) 



where again we note that 9? is arbitrary at r = 0. If the symmetry axis is regular, 7' and ip' 
vanish there, and thus the above quantities vanish. In this case, this curve can be a geodesic, 
since d 2 x/dX 2 and d 2 y/dX 2 vanish, as found from the geodesic equations. By contrast, if a 
spacetime singularity appears on the symmetry axis, 7' and ip' do not vanish there, and 
therefore the above quantities become indefinite, due to the arbitrariness of tp. This implies 
that the curve x = y = = dx/dX = dy/dX cannot be a geodesic at the spacetime singularity, 
since the geodesic equations are not well-defined there. However, it should be noted that 
the geodesies along r = cannot be extended across the spacetime singularity as a causal 
geodesic, but they can be extended as a causal curve, since the metric tensor is defined 
even at the spacetime singular ityj*3 In this sense, this singularity is very weak, like singular 
hypersurfaces, as discussed below. 

5.3. Comparison of the present spacetime singularity with others 

Here it is worthwhile considering the case that r = is conically singular, i.e., e 27 7^ (3 2 
at r = 0. Also in this case, we can uniquely determine a geodesic k + for a given causal 
geodesic fc_ in the same manner as in § 5.2. However, since, in contrast to the cases of 
interest presently, the metric tensor is not defined at the conical singularity, r = 0, we 
cannot conclude the vectors are parallel there. Hence, we cannot claim that the tangent 
vectors to fc_ and k + are parallel to each other at the conical singularity at r = 0, even 
if their components with respect to the Cartesian coordinate system are identical to each 
other in the limit r — > 0. Parallel transport across the conical singularity is impossible, and 
therefore the geodesies that intersect the conical singularity at r = are not extendible as 
continuous geodesies. 

As shown in §§ 5.1 and 5.2, although all of the radial geodesies are complete, the central 
geodesies are not complete in the spacetime of interest. Thus, this spacetime is geodesically 
incomplete. However, this spacetime singularity is not too strong to treat within the frame- 
work of general relativity. The same situations can also be realized in the case of singular 
hypersurfaces, which can be treated within the framework of general relativity with a clear 
physical meaning. 14 ) To make the discussion clear, let us consider a static spherically sym- 
metric shell with infinitesimal thickness but finite mass, which constitutes a singular timelike 
hypersurface. From a physical point of view, such a thin shell is supported by the tangential 

*' The same is true for the null dust case. Two of the present authors have argued that a spacetime 
exhibiting the gravitational collapse of cylindrical null dust is geodesically complete if the passing-through 
condition is imposed, 34 - ) but this argument is not exact. 
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pressure against the gravitational force action to collapse it. We assume that it exists in 
a vacuum spacetime. By the assumed symmetry, the inside of the singular hypersurface is 
described by the Minkowski geometry, 



ds 2 = -dt 2 + dr 2 + r 2 (d6 2 + sin 2 0<kp 2 ), 



(5-17) 



whereas the outside is described by the Schwarzschild geometry, 



-i 



ds 2 = - 



dr 2 + r 2 (dd 2 + sin 2 Odtp 2 ), 



(5-18) 



where r g is a non-vanishing constant smaller than the areal radius of the shell, r = r s h e ii- 
Then, consider curves with constant spatial coordinates r, 9 and ip. They are timelike 
geodesies in the Minkowski domain, but they are not in the Schwarzschild domain. However, 
we cannot determine whether or not a timelike curve attached to this spherical shell, r = 
^sheib is a timelike geodesic, because the connection is not defined on the spherical shell. This 
situation is similar to that of the central geodesies in the present case. 



In order to investigate the strength of the curvature at the spacetime singularity, r = 0, we 
consider timelike geodesies that intersect the spacetime singularity. The curvature strength 
of spacetime singularities was defined in the hope that the weak convergence of geodesic 
congruences would reveal the extendibility of the spacetime in a distributional sense. 33 ^ In 
this context, Tipler defined the strong curvature condition 35 ^ and Krolak defined a weaker 
condition called the limiting focusing condition. 36 ^ 

As mentioned above, non-radial geodesies cannot intersect the spacetime singularity, and 
therefore we consider only radial and central ones. The components of the radial or central 
geodesic tangent e^L are given by 



Then we have an orthonormal frame parallely propagating along this timelike geodesic, 



§6. Curvature strength 




(6-1) 




(6-4) 



(6-3) 



(6-2) 
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where 



(6-5) 



The strength of the tidal force measured by an observer moving along this timelike geodesic 
is given by the components of the Riemann tensor with respect to the above tetrad basis, 
R(o) (A)(0)(B) (A, 5 = 1,2,3). 

We give the coordinate basis components of the Riemann tensor in Appendix D. Using 
those components, it can be seen that the only -R(o)(2)(o)(2) becomes infinite at the spacetime 
singularity, r = 0, and it is given by 



R 



1 



(0)(2)(0)(2) = 
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ip' - 7> - 7^' + ip'ip + ^7' + ^-7 - ^ 
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( dr 
\d^ 



(6-6) 



As in the previous section, the affine parameter r is adjusted so that the geodesies intersect 
the singularity at r = 0. 

Before the singularity formation, f3', 7' and ip' vanish at the symmetry axis, r = 0, and 
thus along the central geodesies, -R(o)(2)(o)(2) is given by 



r — 0-2(7—0) 

«(o)(2)(o)(2) - e 

_ e -2 7 p2 



$ - 7^ - 7V - ^ 2 + 2^ /2 + ^(7 + 2^) + ^ (7' - 2^') - 











(6-7) 



where we have used the normalization of the geodesic tangent. Because the metric functions 
f3, 7 and ip are regular with respect to t and r on the symmetry axis, r = 0, and because 



(3 must be positive, we can easily see that R 



(0)(2)(0)(2) 



is finite in the limit r — > 0— . It 



suddenly diverges just when the conserved densities D and become positive at r = 0. 
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Thus, neither the strong curvature condition nor the limiting focusing condition is satisfied 
along the central geodesies. 33 ^ 

We can see from Eq. (I6-6P that in the case of a radial geodesic with dr/dr < 0, the 
component -R(o)(2)(o)(2) takes the form 



%>)(2)(o)( 2 ) = ^ + J(t), (6-8) 



where / and J are C l functions of the affine parameter r. Thus we have 

I(t) 



Z :=limT|# (0) ( 2)( o)( 2 )| = lim 



dr \ ' 
dr 



(6-9) 



where we have used l'Hopital's theorem in the final equality, 



lim- = limf^) . (6-10) 



t^o r \dr r 

Further, by l'Hopital's theorem, we have 

lim 5 df 5 df \ R W( 2 )(Q)(2)( f )\ _ Um II 2 ^|^(Q)(2)(o)(2)(r)| _ ^ 
r^o r In t — t r^o lnr 

where C\ and C2 are negative constants. Thus we find 

lim I df j df I -R( )(2)(o)(2) I = Z lim(r lnr - r) = 0. (6-12) 

T ^°J Cl Jc 2 

The above equation implies that this spacetime singularity is not the strong curvature sin- 
gularity with respect to the timelike geodesies. 33 ^ We can see that the same is true for null 
geodesies by using the same procedure as above. 

However, this spacetime singularity satisfies the limiting focusing condition, since from 
Eq. (16- 11 1) , we have 

c/f| J R (0 )(2)(o)(2)(^)| = Z lim lnr = oo. (6-13) 



J Co 



'(-'2 

In this sense, this singularity is slightly stronger than singular hypersurfaces. This is an 
example of a spacetime that is extendible across the spacetime singularity which satisfies the 
limiting focusing condition. 33 ) 

§7. Summary and discussion 

We have studied the gravitational collapse of a cylindrical thick shell, in other words, a 
hollow cylinder, composed of dust matter and presented a physically reasonable boundary 
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condition at the resultant string-like naked singularity formed at the symmetry axis. With 
this boundary condition, the trajectories of dust particles can be extended across the naked 
singularity, and thus the dust matter can be regarded as a cold gas of collisionless particles of 
infinitesimal mass; the collapsing dust matter passes through the naked singularity created 
by the dust itself without any interaction between individual constituent particles. When 
the dust matter leaves this region, the naked singularity disappears, and the regularity at 
the symmetry axis is recovered. We performed numerical simulations under this passing- 
through boundary condition and showed that the picture described above is accurate. The 
obtained spacetime should represent a phenomenological description of such a physically 
realistic system of collisionless particles. 

This naked singularity is a scalar polynomial curvature singularity. The strength of the 
curvature at the naked singularity was investigated, and we found that it does not satisfy 
the strong curvature condition defined by Tipler but that it does satisfy the limiting focusing 
condition defined by Krolak. A causal geodesic along the symmetry axis, r = 0, which is 
called the central geodesic, cannot be extended across the spacetime singularity as a causal 
geodesic, but it can be extended as a causal curve. All the other causal geodesies that 
intersect the naked singularity can be extended across the naked singularity. Because the 
central geodesies have zero measure in the space of solutions for the geodesic equations, we 
can conclude that the extended spacetime is complete for almost all of the causal geodesies. 
The present results imply that gravity produced by a cylindrical thick shell of dust matter 
is too weak to bind the shell, even if it engenders a curvature singularity. We conjecture 
that this weakness of gravity is the origin of its nakedness. Such a naked singularity does 
not lead to serious problems, because we can describe this system with a C l ~ spacetime, as 
in the case of singular hypersurfaces. 

In the case of a solid dust cylinder, the spacetime singularity formed by its gravitational 
collapse will be stronger than in the case of a hollow cylinder for the following reason. A 
regular solid dust cylinder should satisfy the relations 

D u 
lim — > 0. and lim — = u c , (7-1) 

r^O r r— >0 r 

where u c is a real number. Because a collapsing solid dust cylinder will have negative u c , u 
will be monotonically decreasing in a sufficiently small neighborhood of the symmetry axis, 
r = 0. This implies that a caustic of dust matter appears at r = within a finite time, and D 
diverges there. The resulting spacetime singularity will be of the third kind defined in § 3.2. 
In this case, provided that the boundary conditions A3- 19)) — (13-21)) are imposed, 7' and ip' 
diverge at r = 0, and therefore the metric tensor cannot be C x ~ for any extended spacetime. 
Thus, we need another prescription to deal with this spacetime singularity. Here, it is worth 
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noting that if an infinitesimal cylindrical portion of dust matter is removed from around 
the center of this solid dust cylinder, it becomes a cylindrical thick shell, and hence the 
spacetime singularity becomes so weak that almost all causal geodesies are complete. This 
means that the strength of the spacetime singularity comes from the infinitesimal central 
portion of the dust cylinder, which seems to be insignificant from a physical point of view. 
Therefore, the spacetime singularity formed by a solid dust cylinder also seems to be so weak 
that we could find an appropriate prescription to treat this spacetime singularity. This issue 
will be discussed elsewhere. 

The Shapiro- Teukolsky singularity, which was reported to be a naked singularity formed 
in the spindle collapse of collisionless particles, might be similar to that studied in the present 
case, except in the neighborhoods of the poles of the spindle matter distribution. In general, 
if the naked singularity formed by the highly elongated gravitational collapse is so weak as 
to allow C l ~~ extension of spacetime, it can also be regarded as phenomenological, employing 
coarse-graining of 'microphysics', such as infinitesimally massive collisionless particles. 

In the present investigation, we did not continue the numerical simulations until nearly 
equilibrium configurations were realized. Due to the pressureless nature of dust matter, shell 
crossing singularities must appear. In the numerical scheme adopted in this paper, the dust 
matter is treated as a system composed of collisionless particles, and thus, in principle, we 
are able to follow the dynamics even after the appearance of shell crossing singularities in 
our numerical code, as long as they do not appear on the symmetry axis, r = 0. This system 
might generate a large amount of gravitational radiation, as in the case of the adiabatic 
gravitational collapse of a cylindrically symmetric ideal gas. 37 -* Investigating this point is 
also a future work. 
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Appendix A 

The Positivity of f3 

In order to demonstrate the relation f3 > 0, we invoke the following theorem proven by 
Hayward: 27 ) If the dominant energy condition is satisfied, the C -energy is non- decreasing 
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in outward achronal direction in untrapped regions. The dominant energy condition for the 
matter of the stress-energy tensor is the condition that T^ V V^ be a causal vector for an 
arbitrary causal vector V^. Here, untrapped means that the relation R 2 — R' 2 < holds, 
where 

R ■= y/9w>9zz = r/3 (A-l) 

is the area element in the 2-space labeled by t and r. The C-energy C is a quasi- local energy 
per unit Killing length along ^-coordinate included within the cylinder of radius r, 38 ) 



C7--1 
C - 8G 



l + e- 2 \R 2 -R' 2 ) . (A-2) 



By the above definition, a non- decreasing C-energy for the outward achronal direction implies 
the relation 

e- 2 \R 2 - R' z ) >0, (A-3) 



where ^ is an arbitrary outward acrhonal vector. 

In the present paper, we consider a spacetime that is identical to the Levi-Civita space- 
time in the spacelike asymptotic region [see Eq. (|4-8|) ]. Since the Levi-Civita spacetime is 
everywhere untrapped, the spacetime of interest has an untrapped spacelike asymptotic re- 
gion. Thus, the inequality (1A-3I) implies that the spacetimes that we consider are everywhere 
untrapped. We have 

R 2 - R' 2 = r 2 p 2 -{(3 + r(3') 2 — ► -f3 2 \ r =o for r — > 0. (A-4) 



Because we have (3\ r =o = e 7 | r= o > 0, the symmetry axis must be untrapped. Equation ( 1A-3I) 
and the above equation lead to 

R' 2 > [3 2 \ r=0 + R 2 . (A-5) 

Since R is a regular function, the above inequality implies that R' > everywhere because 
i?'| r=0 = j3\ r= o > 0. Therefore, R is positive for r > 0, and this means that (3 is positive in 
the physical domain. 

Appendix B 

Non-Radial Geodesies Near the Singularity 



In this appendix, we show that non-radial geodesies cannot intersect the spacetime sin- 
gularity at r = 0. Using Eq. (15-61) . the normalization of the geodesic tangent vector is given 
by 

l 2 < 2 '- 

(It J \dr ) 



e -2(7-V>) 



r 2 (3 2 



+ e - 2 ^P 2 = x , (B-1) 
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where x = ~ 1 f° r timelike geodesies, x = 1 f° r spacelike geodesies and x = for null 
geodesies. Using this normalization condition, we rewrite the r-components of the geodesic 
equations (15-31) in the form 



d 

(It 



dr\ L 2 e 2 ^ f(t 



dr 



r 3 P 



+ 



where 



f(t, r) = L 2 e^ (t - ^ + r V (x ~ e^P 2 ) - rY ( X - 2e^P 2 ) . 



:b-2) 



(B-3) 



Because the metric variables /3 (> 0), 7 and ^ are regular functions of t and r (> 0), there 
exists a positive real number e such that 



sup {rf(t, r) : < r < e} < 1 

holds. Thus, if we consider a positive real number C satisfying 

sup{(L/3e 7 ) 2 : < r < e} 
1 — sup {rf(t, r) : < r < 5} ' 



(B-4) 



C < 



then we have 



d_ 

dr 



,2(7-vo; 



<ir / r 3 



;b-5) 



(B-6) 



Because we consider the future directed causal geodesies approaching the naked singularity 
at r = 0, dr/dr must be negative. Then, multiplying both sides of the above inequality by 

e 2 ^~^dr/dr, we have 



d_ 

d^ 



dr 



_ P 4( 7 -v) 1 ^_ \ _|_ r. 



< 



(B-7) 



for < r < e. From the above equation, there should be a positive constant 7^ such that 

1 



Thus, we have 



2 IdrJ + r 2 



(B-8) 



(B-9) 



The above inequality implies that non-radial geodesies cannot intersect the naked singularity 
at r = 0. 
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Appendix C 

Christoffel Symbols in the Cartesian Coordinate System 



The Christoffel symbols in the Cartesian coordinate system are given by 



r t * = 7 - ip, r} x = (Y - ip') cos ^ r\ y = (y - ip') sin ^, 

= f3(3e~ 21 sin 2 + 7 cos 2 ip - ip (cos 2 </? + /3 2 e~ 21 sin 2 </?) , 



r = 

xy 



/3(3e~ 2 ^ -7 + ^(1- /3 2 e" 27 ) 



smy2cos<^, 



= /3/3e 27 cos 2 + 7 sin 2 <p — ip (sin 2 ip + /5 2 e 27 cos 2 ip) , 

r* = {i-i>') cos ip, 
r £ = | sin 2 up + 7 cos 2 </? - ^, = 



— — 7 | sm ^ cos r \ 







(2 - /3V 27 ) sin 2 up + 7' cos 2 </? - ^ { 1 + (l - V 27 ) sin 2 y?} 



+ -(f-/5V 27 ) sin 2 y9 



COS </?, 



J1X 

xy 



P 



[p 2 e 27 cos 2 ip — cos 2</?) + 7' cos 2 (p — ip' (f3 2 e 27 cos 2 </? + sin 2 ip) 



J-iX 

yy 



■{\ — f3 2 e 27 ) cos 2 (p simp, 
~ [p 2 e- 21 cos 2 up + 2 sin 2 y?) + 7' sin 2 ^ + ip' ((3 2 e~ 2 ~< cos 2 ^ + sin 2 ip) 



H — (l — (3 2 e 27 ) cos 2 p> cos 

r* = -^'e-^+^cosv?, 

^ = (7-^) sin^, 

r te = - ( ^ - 7 J sin </? cos </?, r t y y = ^ cos 2 y? + 7sin 2 up-ip, 



py 



~ (/? 2 e -27 sin 2 up + 2 cos 2 y?) + 7' cos 2 ip + ip' (/3 2 e -27 sin 2 ^ + cos 2 <p) 

i 



+ -(i-/5V 27 ) sin 2 y9 



sm ip, 



py 

xy 



^— {p 2 e 27 sin 2 ip + cos 2ip) + 7' sin 2 ip — ip' {J5 2 e 27 sin 2 ip + cos 2 </?) 
- (l - /3 2 e~ 27 ) sin 2 y? cos </?, 
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py — 

yy 



P (2 - fie-^<) cos 2 <p + i sin 2 <p - {l + (l - /3V 2 ^) cos 2 y?} 



+ i (1 - /3V 2 ^) cos 2 y9 



S1IK/9, 



We can easily see that if 7' and ip' vanish at r = 0, all of the Christoffel symbols are 
single- valued. 



Appendix D 

Riemann Tensor 



The components of the Riemann tensor with respect to the coordinate basis are given by 



rytr 
JX tr 



R tz 



tz 



tz 



R 

R rz 



R 



ipz 



ipZ 



+e 2 ^ (^-Y-^+r), 



+e z 



= — e 



2(^-7) 



^ = 



3 2(V-7) 



(y _ 7 ^ _ 7 y + ^,2 + 2 ^/2 j ; 

^ (7' - + $ - 7^ - 7V + + 1 (7 + v>) 



+ 7j(7-^)-^ 



- _p 2 (V>-7) 



(3 (3' (3' 

- I 7 - ^ I + "0' - 7V - W + "0V + ^7 + - 7j 
+ +^-7^-7y + ^ 2 + ^(7-^) 



with all other components vanishing. 
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